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ABSTRACT 

Context. Several problems regarding the process of galaxy formation are nowadays still open. One of them is the role played by 
Active Galactic Nuclei (AGN) phenomena in contributing to galaxy build-up and in particular to Star Formation (SF) quenching. On 
the other hand, the theory of AGN formation predicts these phenomena to be correlated with the host-galaxy environment, therefore 
opening for links between SF quenching, environment and AGN phenomena in the galaxy formation and evolution paradigm. 

Aims. This work is focused on the study of the correlation between environmental density, radio AGN presence and the probability of 
a galaxy to host a radio Active Galactic Nucleus. 

Methods. Using data from the photometric COSMOS survey and its radio 1.4 GHz follow-up (VLA-COSMOS), a sample of radio 
AGNs has been defined. The environment has been studied throughout the use of the richness distributions inside a parallelepipedon 
with base side of 1 Mpc and height proportional to the photometric redshift precision. Richness distributions have been compared as 
a function of both the redshift and the relative evolution of the stellar masses of galaxies and AGN hosts up to z = 2. 

Results. Radio AGNs are always located in environments that are significantly richer and denser than those around galaxies in which 
radio emission is absent. Therefore, the environment seems to enhance the probability of a galaxy to host a radio AGN. Moreover, a 
distinction between high power and low power radio AGNs leads to the conclusion that the significance in the environmental effect 
is maintained only for low power radio sources. By studying the evolution of stellar masses it is possible to conclude that radio AGN 
presence is a phenomenon related to quiescent galaxies up to z = 2, with a significant increase of the fraction of quiescent galaxies 
hosting a radio AGN with decreasing redshift. Hints of an environmental effect are present as well. 

Conclusions. The results found with this work lead to conclude that denser environments play a significant role in enhancing the 
probability of a galaxy to host a radio AGN and in particular low power ones. 

Key words. Galaxies: active - Galaxies: formation - Galaxies: evolution - Galaxies: mass function - Galaxies: star formation - 
Galaxies: clusters: general 


1. Introduction 

It is well known that the morphology of galaxies is related with 
the local density in clusters (Dressier 1980) Moreover, luminos¬ 
ity functions (Zucca et al. 2009), colours (Cucciati et al. 2010a) 
and stellar masses (Bolzonella et al. 2010) show dependencies 
up to redshift z ~ 1 with the density field. Nevertheless, the 
mechanism that transforms star-forming galaxies into quiescent 
ones in dense environments is not yet well understood. Recently, 
Peng et al. (2010) suggested a model where the star formation 
quenching is separable in a mechanism depending on the galaxy 
mass and on the density; the first mechanism acting principally 
on massive objects at the centre of halos, while the second re¬ 
garding satellite galaxies. 

The mass quenching can be related to the so called “feed¬ 
back”. This term simply indicates an energy input which causes 
an early stop of the star formation in a galaxy. This energy input 
was first introduced to reconcile cosmological simulations with 
observations (see for example Granato et al. 2004; Croton et al. 
2006) and since then it has been naturally associated with the 
activity of an Active Galactic Nucleus (AGN). 

By noting that X-ray and radio AGNs have different clus¬ 
tering properties, Hickox et al. (2009) suggested that the cen¬ 
tral AGN co-evolves with the host galaxy: while the host galaxy 
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transforms from a star-forming to a quiescent one, the AGN 
passes from a quasar. X-ray emitter phase to a radio-galaxy 
one. These transformations happen at earlier epochs for halos of 
higher mass. This evolution of AGN type with galaxy transfor¬ 
mation could be the reason behind the estimate that radio AGNs 
constitute ~ 10% of the whole Active Galactic Nuclei population 
(see for example Hickox et al. 2009), although projection effects 
and the limiting flux of current radio surveys may also play a 
role. The motivation of such a small fraction of radio to total 
AGNs is still not fully understood. Since the analysis conducted 
in this work focuses only on radio AGN samples, in further text 
it will be explicitly referred to this kind of sources. 

In the scenario developed by Hickox et al. (2009) and Peng 
et al. (2010), the environment is still related to the star forma¬ 
tion quenching of massive galaxies in an indirect way, in the 
sense that massive galaxies (where the mass quenching is effec¬ 
tive) reside preferentially in high density environments, where at 
low redshifts early-type galaxies dominate (Quadri et al. 2012; 
Chuter et al. 2011) Moreover, it was already known that many 
radio AGNs reside in early-type galaxies (Ledlow & Owen 
1996) and that the probability for a galaxy to host an AGN is 
increasing with stellar mass (Bardelli et al. 2009) Also, Bardelli 
et al. (2010), on the basis of the zCOSMOS 10k sample (Lilly 
et al. 2007, 2009) showed that the fraction of radio active early- 
type galaxies is an increasing function of local density. 
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However, what it is not clear is what is meant when envi¬ 
ronment is considered, since its definition is related to the spa¬ 
tial scale studied. In fact, environments estimated with nearest 
neighbour methods consider different scales at different redshifts 
in a non trivial way. For example, at least at low redshifts. Best 
(2004) claims that the larger scale (i.e. the scale typical of groups 
and clusters) is more important than the smaller scale (typical of 
galaxy pairs or companions) to determine the AGN radio emis¬ 
sion. 

In this work the environment of radio sources of the VLA- 
COSMOS survey (Schinnerer et al. 2007), cross-identified with 
the COSMOS photometric redshift sample (Ilbert et al. 2009), 
is explored. With respect to the Bardelli et al. (2009) work, this 
sample allows an increase in both redshift range and statistics, 
but at the price of having a larger error in the redshifts (photo¬ 
metric versus spectroscopic redshifts), which excludes the possi¬ 
bility of a small spatial scale definition for the environment. The 
scale of the environment is fixed at 1 Mpc, in order not to depend 
on the scale-redshift relation. 

The data and the method used to define the radio AGN sam¬ 
ple are described in sections [^andj^respectively. In section|^the 
method to compare the overdensity richness distributions is pre¬ 
sented and in sections|^and|^the results for various sub-samples 
are discussed. In section [7] a correlation with known clusters and 
groups catalogues is attempted and in section the integrated 
mass functions of radio AGNs and normal galaxies are presented 
and analysed. Conclusions are drawn in section Throughout 
the text, masses are expressed in units of Solar masses, while 
radio luminosities are expressed in SI units. The assumed cos¬ 
mology is flyi = 0.7, flm = 0.3 and Hq = 70 km s“^ Mpc“^ For 
the SED fitting procedure, in order to derive stellar masses and 
Star Formation Rates, Ilbert et al. (2010) used SED templates 
produced with a Chabrier (2003) IMF. 


2. Data 

In this work, two main datasets have been used. Optical data 
(with photometry in 30 bands from UV to mid-IR) have been 
used to estimate photometric redshifts for all the sources, to de¬ 
termine the environment around every object in every sample 
and to extract Control samples as explained later in the text. 

Radio data at 1.4 GHz, instead, have been used to create the 
radio source and AGN samples, around which the environment 
has been studied. 


2.1. Optical data 

The data at optical wavelengths are constituted of the photomet¬ 
ric sample of Capak et al. (2007), used in the construction of the 
version 1.8 of the photometric redshift catalogue of Ilbert et al. 
(2009) This sample is a compilation of photometric data taken 
from the COSMOS survey (Scoville et al. 2007), which covers 
a field of an area of about 1.4“ x 1.4“ and is centred at R.A. 
(J2000.0) = 10^00"^28.6" and Dec (J2000.0) = -r02° 12'21.0". 

In Ilbert et al. (2009), photometric redshifts have been calcu¬ 
lated through a Spectral Energy Distribution (SED) fitting pro¬ 
cedure, using 30 broad-, intermediate- and narrow-band filters 
from UV to mid-IR frequencies (see their Table 1). The fit has 
been performed with slx^ minimization algorithm on a template 
library using the Le Phare code (S. Amouts & O. Ilbert). 

The redshifts obtained with the described procedure have 
been compared with those obtained using the zCOSMOS spec¬ 
troscopic survey (Lilly et al. 2007) 
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The accuracy of the photometric redshifts {zp) is estimated 
to be 

o'^z|{\+z) = 0-06 ( 1 ) 

This value will be assumed throughout this work. The error 
has been conservatively selected as the maximum value of those 
obtained as a function of both magnitude and redshift in Ilbert 
et al. (2009, see their Figure 9), taking into account the depth of 
the optical catalogue used for the analysis. In order to check the 
consistency of such a choice, the effect of degrading the value to 
<^Az/(i+z) = 0.1 ^ 0.2 for high redshift sources has been tested. As 
explained further in the text, it has been found that it does not 
change significantly the results of this work. 

Other sources of uncertainties in the photometric redshifts 
determination which must be kept in mind are the so called 
catastrophic errors. The definition of catastrophic errors is ap¬ 
plied to those sources for which the redshift calculation fails in 
the form of 

> 0.15 (2) 

where is the spectroscopic redshift. It is estimated that the 
fraction p of catastrophic errors rises from 1% to 20% going 
from C < 22.5 to C ~ 24. 

Only objects that are not in masked areas and with magni¬ 
tude C < 26.5 are considered. This magnitude limit approxi¬ 
mately corresponds to the point where the magnitude-number 
counts histogram stops increasing. At this limit the average num¬ 
ber of filters used for the photometric redshifts determination is 
18. Although the effect of such a deep cut in limiting magnitude 
is to increase the uncertainties in the physical quantities derived 
from the optical photometry, it has also the advantage of greatly 
increasing the statistics for the environment estimate. For this 
reason, this cut will be assumed for all samples in this work. As 
explained extensively in section]^ tests have been performed by 
introducing brighter magnitude cuts, from C < 25.5 to C < 24, 
to check the problematics related to including in the analysis also 
very faint objects. These tests showed that the magnitude cut has 
no major effects on the results found with this work. 

From the optical catalogue a sample of 823 939 sources with 
optical data and measured photometric redshift has been ex¬ 
tracted (hereafter “O sample”). Star Formation Rates (SFR) and 
stellar masses (M*) are both derived through SED fitting with 
population synthesis models (see Ilbert et al. 2010), together 
with other physical parameters. 

2.2. Radio Data 

The radio data are taken from the VLA-COSMOS Large Project 
version 2.0 catalogue (Schinnerer et al. 2007, 2010), whose ob¬ 
servations have been carried out with the Very Large Array in its 
A configuration. This grants a resolution of about 1.5" at Full 
Width Half Maximum (FWHM) in the observation wavelength 
of 1.4 GHz. A total of 23 pointings was necessary to cover the 
full COSMOS field, for 240 hours of observation, performed be¬ 
tween 2004 and 2005. The observations reached a sensitivity of 
about 11 yuJy r.m.s. and the 5cr detection limit catalogue contains 
2417 radio sources, 78 of which have been classified as multi- 
component radio sources. 

The search for the optical counterparts has been done us¬ 
ing the photometric redshift sample (version 1.8) by Ilbert et al. 
(2009), without limitation on the optical magnitude, while for 
the radio catalogue the VLA-COSMOS Large Project version 
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Table 1. Size of every sample and redshift range of every bin used in 
the data analysis. 


Bin 

z 

R sample 

MR sample 

AGN sample 

Bin 1 

o 

o 

532 

425 

119 

Bin 2 

0.7-1 

320 

290 

100 

Bin 3 

1-2 

575 

503 

53 



z 


Fig. 1. Redshift Distributions. The black solid line refers to the O sam¬ 
ple, the red dashed one to the R sample, the cyan dot-dashed line refers 
to the MR sample and the green triple dot-dashed one refers to the AGN 
sample. 


2.0 has been used, but without the 78 multi-component radio 
sources (2339 radio sources). 

The method used for the optical identification is the likeli¬ 
hood ratio method, described in Sutherland & Saunders (1992), 
Ciliegi et al. (2005) and Brusa et al. (2007) For the likelihood ra¬ 
tio analysis it has been adopted a Q value (the probability that the 
optical counterpart of the radio sources is brighter than the mag¬ 
nitude limit of the optical catalogue) of 0.9 and a likelihood ratio 
cut-off of 0.1 (1 - Q). Using a likelihood ratio cut-off equal to 
l-Q ensures that ah the optical counterparts of the radio sources 
with only one identification and a likelihood ratio greater than 
the cut-off value have a reliability greater than 0.5. With these 
thresholds 2069 radio sources with an optical counterpart have 
been found, ~ 26 of which (~ 1.3%) could be spurious positional 
coincidences. Applying to the 2069 identified radio sources the 
same cuts used to define the O sample (objects not in masked 
areas and with magnitude t < 26.5) a sample of 1806 optically 
identified radio souces has been created (hereafter “R sample”). 
The data analysis has, nevertheless, been performed only out to 
z < 2.0, i.e. on a selection from the R sample of 1427 objects. 
The sizes of the various samples can be found in Table while 
a summary of the sample definitions can be found in Table 

In Figure the redshift distributions for the O, R, MR and 
AGN samples (the last two will be introduced further on in the 
text) are compared. It can be seen that the distributions referring 
to radio galaxies peak at lower redshift, while the distribution for 
the optical galaxies is wider in a more extended redshift range. 

The greatest difference in the samples is highlighted in Fig¬ 
ure 1^ which shows the stellar mass distributions for the same 
samples as in the previous plot. It can be easily seen that the 
galaxies showing radio emission (R, MR and AGN samples) are 
found almost only in high mass galaxies (the peak is around 



Log(M*/M_Sun) 


Fig. 2. Stellar Mass Distributions. The black solid line refers to the O 
sample, the red dashed one to the R sample, the cyan dot-dashed line 
refers to the MR sample and the green triple dot-dashed one refers to 
the AGN sample. 


log(M*) ~ 11), while objects from the O sample span a much 
more wide range of masses, with a peak at log(M*) ~ 9. 

This could be due both to an evolutionary effect, that sees 
the onset of radio emission preferentahy in high mass galaxies, 
or to the limiting flux of the VLA-COSMOS which could limit 
the kind of galaxies detected to only the most massive and lu¬ 
minous ones. Aside from the physical process responsible for 
the shape of the stellar mass distributions, the fact that the radio 
galaxy samples are located at higher masses with respect to nor¬ 
mal galaxies must be kept into account, in order not to introduce 
a bias in the conclusions drawn using these samples. For this rea¬ 
son a cut in stellar mass has been introduced, selecting only the 
most massive radio sources. Considering only radio sources at 
log(M*) > 10 the sample size becomes of 1448 objects, (“MR 
sample”). Again, the data analysis has been conducted only on 
a selection of 1218 with z < 2.0. The cut to log(M*) > 10 also 
accounts for the fact that the radio sources without optical identi¬ 
fication could be faint galaxies with low stellar mass. Therefore 
selecting only high mass radio sources should avoid the insur- 
gence of any bias in the results from missing optical identifica¬ 
tions of radio sources. 


3. The AGN sample 

It is well known that radio emission could be due, in principle, 
both to Star Formation (SF) and to AGN phenomena. Several 
methods have been suggested for dividing the two populations 
where each of the two phenomena is dominant. There are several 
recipes for segregating the radio population into sources with 
emission from an AGN or from star formation (see e.g. Smolcic 
et al. 2008; Bardehi et al. 2009; Best et al. 2005, and references 
therein). In this work two ways which share a common hypoth¬ 
esis have been adopted. 

The first way of defining a sample of radio AGNs is through 
a cut in Specific Star Formation Rate (SSFR). For each ob¬ 
ject the SSFR is defined as the ratio of the SFR and the stel¬ 
lar mass, with both quantities derived from the SED fitting pro¬ 
cedure. Quiescent galaxies are selected and it is assumed that 
for this kind of sources ah the radio flux comes from AGN 
activity. For this purpose, the Specific Star Formation Rate 
Stellar Mass (SSFR - M*) plane is analysed (Figure |^. From 
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Table 2. Samples definition. 


Sample Name 
O 
R 

MR 

AGN 

RO 

MO 

QO 


Definition 

Optical objects, t <26.5 
Optically identified radio objects, t <26.5 
Optically identified radio objects, t < 26.5, log(M*) > 10 
AGNs, t < 26.5, log(M*) > 10, log(^) < -11 

yr 

Control sample of R, extracted from O 
Control sample of MR, extracted from O, same M* distribution 
Control sample of AGN, extracted from O, same M* distribution, same SSFR range 



^ -10 
i 

F -12 


-14 


-16 

-18 


+ R Somple 

A AGN Sample 

□ SFR Componson Somple 


iHi «ii:T:Knii|k 


7 8 9 10 11 12 

Log(M*/M_Sun) 



Log(SFR_SED/M_Sun*yr-1) 


Fig. 3. SSFR - M* plane. The black crosses refer to the R sample, the 
red triangles to the AGN sample defined through the SSFR cut and the 
cyan squares to the AGN sample defined through SFR comparison. The 
green lines correspond to the cuts in Stellar Mass and Specific Star For¬ 
mation Rate described in text. In yellow, the O sample is reported for 
comparison. 


Fig. 4. AGN sample definition. The black crosses represent the R sam¬ 
ple, the red triangles represent the AGNs selected through the cut in 
Specific Star Formation Rate, while the cyan squares represent the AGN 
sample defined through SFR comparison. The dashed green line repre¬ 
sents the linear fit relation described in text, while the solid blue line 
shows the S FRradio - 5 FRs ed relation. 


this plot, it is clearly visible the locus of star-forming galax¬ 
ies as an horizontal region within a narrow range of SSFR 
(-11 < log(SSFR/yr~^) < -1). On the right of the figure, a 
vertical band located at 10 < log(M*) <12 represents the galax¬ 
ies which are switching off their star formation (“the dead line”). 
As AGN sample it has been chosen to select radio sources with 
log(M*) > 10 and \og{SSFRIyr~^) < -11. With these cuts, 
a sample of 276 radio Active Galactic Nuclei has been defined 
(272 at z < 2). 

Although efficient, this method adopts an extreme hypoth¬ 
esis, since it assumes that all radio flux is coming from radio 
AGN activity in quiescent galaxies. The cuts applied to the sam¬ 
ple may actually result inaccurate as more reliable methods to 
discriminate between radio activity from AGN phenomena and 
from star formation can be developed. 

One of these is described in Bardelli et al. (2010) It con¬ 
sists in the comparison between Star Formation Rates as resulted 
from the SED fitting of Ilbert et al. (2010) and as estimated from 
the radio power as done in Bell (2003) The hypothesis is that all 
radio flux exceeding that given by the SFRsed is due to AGN 
activity. 

The S FRRadio has been computed using the Bell (2003) for¬ 
mula (which is derived for a Salpeter IMF), with a correction 
of ~ 2 in the normalization introduced to be consistent with 
Herschel-PEP data (Bardelli et al. in prep.). This normaliza¬ 
tion correction also takes into account the difference between 
IMFs used in the SFR derivation from optical SED fitting and 
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in the formula from Bell (2003) Note that in this particular 
work the normalization is less important, since only the iden¬ 
tification of the locus of star forming galaxies matters. The com¬ 
parison between the two star formation rates is reported in Fig¬ 
ure In this figure it is clearly visible an elongated region due 
to the star-forming galaxies with a slight zero point shift and 
tilt. In order to fit this relation, only those objects that have 
log(5' FRradio) ^ log(5' FRs ed)f2 have been considered. The lin¬ 
ear fit resulted in log(5'= 1.176 xlog(5'FR5£/))-F 0.685. 
The distribution of the distances from the best-fit line has then 
been fitted with a Gaussian distribution, obtaining a dispersion 
of 1.40. In order to be very conservative, all the objects that lie at 
a distance greater than 3(T from the fit line are defined as AGNs. 
With this method, 164 sources are selected (154 at z < 2). 

It has to be noted that this second method is more conserva¬ 
tive, since the assumptions are less extreme: it is now possible 
for every galaxy, even for the quiescent ones, that part of the 
radio emission could be due to residual star formation. The dis¬ 
crimination is set to decide which of the two phenomena (radio 
AGN activity or the formation of stars) prevails on the other. 
With this method it is increased the probability of correctly se¬ 
lecting a radio AGN source (purity), in spite of completeness, 
since some radio AGN sources may not be individuated by the 
selection criteria. Figures and both show that the radio AGN 
samples derived with the two methods overlap. 

Since the infrared emission is a better tracer for the SFR, 
a more rigorous approach would be to apply the method de- 
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scribed above using the SFR derived from IR data, rather than 
from optical ones. For this reason the R sample has been also 
matched to the DRl catalogue of the Herschel PACS Evolution¬ 
ary Prob^(Herschel-PEP). A description of the survey, obser¬ 
vational strategies and data reduction techniques may be found 
in Lutz et al. (2011), while SFR and physical quantities have 
been obtained through a SED fitting procedure performed with a 
modified version of the MAGPHYS code (see Berta et al. 2013; 
Delvecchio et al. 2014, and references therein for a detailed de¬ 
scription) A total of 923 sources have been found matching the 
R sample, 45 of which are in the AGN sample and are all situated 
at lower SFR values compared with other Herschel sources. 

From the comparison of the S FRRadio with the SFR from the 
IR data, derived through a procedure of SED fitting that allows 
for the correction of the AGN contribution to galaxy emission 
{S FRjr), a sample of 83 AGNs is created. The method used is 
equal to the one applied in previous paragraphs: the S FRRadio 
and S FRjr are compared and the bulk of the relation is linearly 
fitted in order to isolate sources with an excess of S FRRadio with 
respect to S FRjr. The distribution of the distances of each object 
from the best-fit line is fitted with a Gaussian distribution and 
only sources residing farther from the best-fit line than 3(T are 
defined as AGNs. 

The environmental analysis has been performed on these 
sources in the same way as for the other AGN samples. The envi¬ 
ronmental segregation effect is recovered, but unfortunately the 
sample is too small to grant significativity to the results. 

Extensive tests have been performed to check whether the 
various methods of defining AGNs yield differences in the re¬ 
sults found in this work. No significant discrepancies emerge 
from the tests. The AGN sample defined through the cuts in 
SSER and stellar mass will be used since it permits a slightly 
larger statistics, but all the conclusions that will be derived in 
the following remain valid also if the radio AGNs are selected 
through the S FRRadio and S FRs ed (or S FRjr) comparison. 


4. Method 

The aim of this work is to investigate the presence and proper¬ 
ties of galaxy overdensities in the COSMOS Survey field, espe¬ 
cially at high redshift. Eirst of all, a definition of overdensity is 
needed, since many different ways of defining this quantity can 
be found in the literature (see Scoville et al. 2013; Kovac et al. 
2010; Cucciati et al. 2010b) and new ways are still undergoing 
developement (see e.g. Castignani et al. 2014a) In the present 
work it has been chosen to have the most simple approach by 
counting objects in a parallelepipedon with a base side of 1 Mpc 
(comoving), centred on the considered source, and height twice 
Az = 3 X (TazHi+z) X (1 + Zp). It is useful to define three kinds of 
samples, to which it will be referred in the following. 

Starting sample. This kind of samples is composed of the ob¬ 
jects around which the environment is to be estimated and 
that will be placed at the centre of each parallelepipedon. 
Starting samples will be the R, MR and AGN samples. 

Target sample. These objects are those that are counted inside 
each parallelepipedon and that are used to estimate the envi¬ 
ronment around every galaxy of each Starting sample. The 
O sample will be used as Target sample. 

Control sample. The absolute values and richness distributions 
of overdensities do not have any direct meaning if not com¬ 
pared with those obtained with a Control sample. In order 
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to create it, for each galaxy of a Starting sample the redshift 
has been maintained, but coordinates have been randomati- 
cally extracted. The way in which this is achieved is different 
for every sample and it will be thoroughly explained in sec¬ 
tion [4[T] The environment has then been re-estimated for ev¬ 
ery Control sample using the same Target sample as for the 
corresponding real ones. The richness distributions between 
every sample and its corresponding Control one have been 
compared with a Kolmogorov-Smirnov test. In order not to 
be influenced by only one extraction this process has been 
repeated 100 times. 


The environmental analysis has been performed in three dif¬ 
ferent redshift bins, namely z e [0.0 - 0.7 [, [0.7 - 1.0[ and 
[1.0 - 2.0]. It is important to mention that all the environmen¬ 
tal analysis has been performed both on a 1 Mpc and on a 2 Mpc 
scale. It has been found that all the results exposed in the fol¬ 
lowing sections have their significance greatly increased in the 1 
Mpc compared to the 2 Mpc case. This could be due to the fact 
that when adopting a smaller scale only the central regions of 
the bound structures are considered, especially at high redshfit, 
where clusters and groups were in an earlier stage of formation. 
It is known that in such regions {i.e. at the centre of overdensi¬ 
ties) galaxy formation happens on shorter timescales. Therefore 
any existing relation between processes involved in galaxy evo¬ 
lution and environment (such as the presence of radio AGNs) 
is likely to be better evidenced on a 1 Mpc scale. Larger scales 
(such as the 2 Mpc case) may have the signal diluted by the in¬ 
clusion of regions in structures that are not yet well formed and 
by an increase in the number of interlopers due to the photomet¬ 
ric redshfit error. This is confirmed by the fact that the results of 
this work are much more significant when using the 1 Mpc scale 
compared to the 2 Mpc case. Eor this reason, in the following, 
only the results obtained with an estimate of the environment on 
a scale of 1 Mpc will be discussed. 

As already mentioned in previous sections, several tests have 
been performed to check whether the choice of the photomet¬ 
ric error and of the magnitude cut performed on the sample do 
affect the results of this work. In particular, the environmental 
analysis has been performed by using also a cut on the O sam¬ 
ple at brighter magnitudes, namely P < 25.5, P < 25 and 
P < 24. Only the magnitude cut on the O sample has been 
changed, since no AGN sources have 25 < P < 26.5. In this 
way, while the magnitude distribution for the AGN sample is 
maintained, that of the Control samples and of the sources used 
to trace the environment has been changed. In the first two mag¬ 
nitude cuts the assumption of the photometric redshift error has 
also been degraded to o-azHI+z) = 0.1 0.2 in the third redshift 

bin (z € [1.0 - 2.0]). The significativity of the results is main¬ 
tained for all the combinations, except that for the AGNs in the 
highest redshift bin when using a magnitude cut for the O sample 
of P < 24. This is not unexpected, though, since tracing the en¬ 
vironment using only the most luminous sources and comparing 
AGNs with luminous control galaxies has the effect of compar¬ 
ing objects residing (especially at high redshift) in similar envi¬ 
ronments, therefore reducing the significativity of the signal of 
environmental segregation. It is therefore possible to conclude 
that neither the magnitude cut used, nor the value of the pho¬ 
tometric redshift error assumed have a systematic effect on the 
results of this work. 
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4.1. Control samples definition 

As briefly explained above, the Control samples for the R, MR 
and AGN subsamples have been created with optical data. For 
every galaxy in the R, MR and AGN samples, coordinates have 
been randomly extracted from the O sample in the same red- 
shift bin for every Starting sample. In addition, the following 
prescriptions have been used separately for every subsample. 

R Sample. The Control sample RO has been extracted from O, 
so to have the same number of sources as the R sample in 
every redshift bin. 

MR Sample. The Control sample MO has been extracted from 
O, so to have the same number of sources as the MR sample 
in every redshift bin. Moreover, every galaxy of the Control 
sample has been extracted among those having stellar mass 
in an interval of +0.5 dex from the mass of each galaxy com¬ 
posing the MR sample. In this way, the MO sample has been 
extracted having the same stellar mass distribution as the MR 
one. 

AGN Sample. The Control sample QO has been extracted from 
O, so to have the same number of sources as the AGN sam¬ 
ple in every redshift bin. It has also been extracted having the 
same stellar mass distribution as the AGN sample, with the 
same method used to extract the MO Control sample. Also, 
the galaxies used as a pool for the QO sample extraction 
have been selected to be quiescent (log(5 'S FRIijr~^) < -11), 
therefore the QO sample is located in the same lower right 
region of Figureas the AGN sample. 

5. Results 

In this section the results of the analysis of the environment per¬ 
formed on the various Starting samples are analysed. As ex¬ 
plained before, in order to better evidentiate the differences in 
the environment around radio sources and around galaxies with 
no radio emission (represented by the extracted Control sam¬ 
ples), the overdensity richness distributions have been compared 
through the use of a Kolmogorov-Smirnov test, repeated with 
100 independent extractions of the various Control samples. 
When applied, the KS test gives the probability that two distribu¬ 
tions come from the same parent population. In order to quantify 
the difference in the environments recovered, only richness dis¬ 
tributions that have a KS probability value lower than 0.05 of 
being extracted from the same parent population will be consid¬ 
ered as significative. In Table the number of Control sample 
extractions (out of 100) that are below the significativity thresh¬ 
old of 0.05 are reported, together with the median values of the 
KS probability value distributions (in brackets). 

5.1. R Sample 

At first, the environment is estimated around sources from Start¬ 
ing samples R and RO. It has to be kept in mind that the R sample 
is the whole radio catalogue. This sample does not discriminate 
whether the radio emission comes from star formation or AGN 
activity. The Target sample used is represented by the whole O 
sample. Results of the KS test between R and the 100 RO extrac¬ 
tions are listed in the first line of Table [S] From these values it 
can be seen, already, that the environment around radio sources 
is significantly different from the environment around sources 
with no sign of radio emission, in every redshift bin. The total 
of the Control sample extractions is always below the significa¬ 
tivity threshold, with very low median values. The way in which 
these environments differ is exemplified by Figure 
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In this figure, the overdensity richness distributions for the R 
sample and for all the 100 extractions of the RO Control sam¬ 
ple are shown. The richness distribution of the Control sample 
extraction that corresponds to the median of the KS probability 
value distribution is in evidence. It can be seen that in every red¬ 
shift bin the radio sources are sistematically distributed at higher 
overdensity richness than normal galaxies. 

Although important, this result could be due in principle to 
other effects, such as the different mass distributions or evolu¬ 
tionary status of R and RO galaxies. A deeper investigation of 
the phenomenon is thus needed. 

5.2. M Sample 

It has been shown (see Figure ^ that optical galaxies and ra¬ 
dio sources have different mass distributions. In particular, radio 
sources are preferentially found in high mass galaxies, which in 
turn are known to reside preferentally in high density environ¬ 
ments. This fact could be the actual reason behind the difference 
in environmental density found between sources from the R and 
RO samples, described in the previous section. 

In order to rule out this effect, the sample MR and its corre¬ 
sponding Control sample MO have been created. These samples 
are designed to include only massive galaxies (log(M*) >10) 
and to have the same stellar mass distribution. In this way, being 
the MR sample composed of radio sources and its correspond¬ 
ing Control sample MO of optical galaxies, any difference in the 
environment between the two can be ascribed to the actual pres¬ 
ence of radio emission. 

The second line of Table [^confirms the trend already found 
for the R and RO sample. The environment around radio sources 
and the environment around normal galaxies are significantly 
different, as demonstrated by the fact that all 100 extractions of 
the MO Control sample are below the significativity threshold, 
with very low median values. Again a visual inspection of the 
overdensity richness distributions, as those reported in Figure 
for the MR sample and all of the 100 extractions of the MO 
Control sample, confirms that overdensities found around mas¬ 
sive radio sources are distributed at higher richness values than 
those around massive galaxies without radio emission. 

It is therefore possible to exclude a stellar mass effect behind 
the differences in the environment around normal galaxies and 
radio-emitting ones. Instead these differences are to be ascribed 
to the physical process behind the radio emission and its corre¬ 
lation with environmental density. 

6. The environment of radio AGNs 

The main focus of this paper is the study of the environment 
around radio AGN sources in the VLA-COSMOS Survey field. 
This task has been performed in a similar way also in other works 
(for example Best 2000; Hatch et al. 2014) In this work a general 
approach is adopted and the results are in good agreement with 
the literature. 

The radio AGN sample for environment definition has been 
extracted from the whole population of radio sources in the way 
described in sectionBy using a cut in Specific Star Formation 
Rate, only sources whose radio emission is due to AGN activity 
have been selected and their environment has been determined 
using the O sample as Target sample. 

Nevertheless, the method used to isolate radio AGNs has the 
effect of selecting also radio sources hosted by massive and qui¬ 
escent galaxies. These two kind of objects are both naturally 
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Table 3. Number of extractions of the various Control samples that are below the significativity threshold of 0.05. In brackets, the median value of 
the KS probability value distributions is given. The probability resulting from the KS test is the probability of the two distributions analysed being 
extracted from the same parent population. 


Sample 

0 < z < 0.7 

0.7 < z < 1 

l<z<2 

R vs RO 

100 (6.6 X 10-«) 

100 (7.5 X 10-^) 

100 (6.5 X 10-') 

MR MO 

100 (6.5 X 10“^) 

100 (6.9 X 10-^) 

100 (9.2 X 10“®) 

AGN vs QO 

100 (8.6 X 10-5) 

100 (6.0 X 10-®) 

87 (0.006) 


found in high density environments, therefore leading to a pos¬ 
sible bias in the conclusions drawn. In fact, any difference in the 
environment around radio AGNs compared to optical galaxies 
could be actually due to quiescent galaxies being found pref¬ 
erentially inside overdensities or to massive galaxies being the 
more clustered. In order not to be affected by these problems, 
the QO Control sample has been selected carefully to avoid this 
bias. Galaxies in the QO sample have been extracted in the same 
SSFR range and with the same stellar mass distribution of the 
AGN sample. In this way, every difference in the environmental 
properties can be ascribed to the presence of radio AGN phe¬ 
nomena. 

The results of the comparison of the overdensity richness dis¬ 
tributions are shown in the third line of Table O It can be seen 
that the environment around radio AGN sources is significantly 
different from the environment around normal galaxies, with the 
vast majority of the extractions of the QO sample being under 
the significativity threshold. Also the median values confirm this 
trend, being always under 5%. This is the most important result 
of this work, showing a real correlation between environment 
and the presence of radio AGN phenomena. 

Figure [7] shows the overdensity richness distributions for the 
AGN sample and for the 100 extractions of the QO Control sam¬ 
ple. It can be clearly seen that overdensities around radio AGNs 
are distributed at higher richness values than those around nor¬ 
mal galaxies. It is therefore possible to conclude that the environ¬ 
ment around radio AGNs is significantly richer than around nor¬ 
mal galaxies, with this effect being due to radio AGN presence, 
rather than to radio AGNs being found preferentially in massive 
or quiescent galaxies (see for example Hickox et al. 2009 and 
Ledlow & Owen 1996 ). In other words, galaxies in denser envi¬ 
ronments show an enhanced probability of hosting a radio AGN 
than galaxies in less dense environments. 

As a check on the results it has also been estimated the envi¬ 
ronment around a sample of radio-emitting starforming galaxies, 
selected in the upper right quadrant of Figure These sources, 
characterized by log(M*) > 10 and \og{SSFRIyr~^) > -11, 
have been compared with a Control sample extracted in the same 
SSFR region and with the same stellar mass distribution. An en¬ 
vironmental effect has been found in this sample too, at all red- 
shifts, with star-forming radio sources being in higher density 
environments. The median values of the KS probability value 
distribution, though, are higher by an order of magnitude with 
respect to those in the last line of Tableexcept that in the high¬ 
est redshift bin, even if the sample of star-forming radio sources 
is ~ 10 times greater than the AGN one (450 SF sources against 
53 AGNs). Although the great inequality in statistics between 
radio AGNs and star-forming sources surely plays a major role 
in determining the significativity of the environmental eflTect for 
the latter kind of objects, cause of this are also probably residual 
differences in the Star Formation Rate between the star-forming 
radio sources and the Control sample ones. This is evidenced by 
the fact that, although they are situated in the same SFR range, 
the star-forming radio sources are detected by the VLA instru- 


Table 4. Number of galaxies in each sample for every redshift bin. High 
power and low power distinction. 


Sample 

0 < z < 0.7 

0.7 < z < 1 

1 <z<2 

AGNH and QOH 

6 

16 

19 

AGNL and QOL 

16 

30 

17 


ment while the Control sample galaxies are not. Therefore the 
SFR is higher in the radio-detected sample compared to the Con¬ 
trol one, or there could still be some residual contamination from 
AGN emission, which induces the environmental effect on these 
sources. 


6.1. High power and low power radio-sources 

It has been suggested (see Miley & De Breuck 2008; Chiaberge 
et al. 2010) that low power radio AGNs are often found in high 
density regions (see also Castignani et al. 2014b; Bardelli et al. 
2009; Donoso et al. 2010) In fact, low power radio AGNs are 
used to detect protoclusters at high redshifts (see e.g. Cooke 
et al. 2012; Wylezalek et al. 2013) In the next it is then ex¬ 
plored if the environmental density is actually depending on the 
radio power. To do this, every sample has been divided in a high 
power subsample (log(Li 4GHz) ^ 24.5) and a low power one 
(24 < log(LiAGHz) < 24.5). The reason for the lower limit at 
log(Li 4GHz) = 24 is to avoid the effect induced by the flux limit 
of the VLA-COSMOS Survey. Such an effect is visible in Figure 

m 

This figure shows the dependence from redshift of the ra¬ 
dio luminosity for the R sample. It can be seen that at high red- 
shifts only the most luminous sources will be detected, due to 
the survey limiting flux (effect known as the Malmquist Bias). 
This could lead to a bias in the results based on the sample, since 
more luminous objects could represent a different population of 
radio sources with respect to less luminous ones. This problem 
can be solved by considering only the sources more luminous 
than \og{Li 4GHz) = 24. In this way, a complete sample of radio 
sources up to z = 2 can be created. 

It is worth noting that the limit at \og{Li 4GHz) - 24.5, set to 
divide the high power and low power radio subsamples, roughly 
corresponds to the canonical division between FRI and FRII ob¬ 
jects (Fanaroff & Riley 1974) 

Therefore the samples AGNH and QOH (where H in the 
sample name stands for high power) and the samples AGNL and 
QOL (where L stands for low power) have been created. The 
number of galaxies in each sample is shown in Table 

The results of the analysis of the environment in relation with 
the distinction in radio power are now discussed. In particular, in 
Table the number of Control sample extractions below the sig¬ 
nificativity level of 0.05 are reported, together with the median 
values of the KS probability value distribution. 

It is possible to see that the high power and low power sam¬ 
ples have different KS probabilities: the environmental effect re- 
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Table 5. Number of extractions of the various Control samples that are 
below the significativity threshold of 0.05. In brackets, the median value 
of the KS probability value distributions is given. Low power and high 
power radio source distinction. 


Sample 

0 < z < 0.7 

0.7 < z < 1 

l<z<2 

AGNH vv QOH 
AGNL vv QOL 

9 (0.08) 

93 (0.002) 

43 (0.07) 

99 (0.0006) 

17 (0.25) 
37 (0.08) 


mains present only for low power radio sources. In fact, in Table 
[^it can be seen that the number of extractions below the signi¬ 
ficativity level is fairly low for the high power radio AGNs, never 
reaching even 50% of the total, and the median KS probability 
value is always greater than 0.05. On the other hand, for low 
power radio AGNs the number of extractions below the signi¬ 
ficativity threshold is close to 100%. Only in the highest redshift 
bin the value is quite low (37%), but it is nevertheless larger than 
the corresponding value in the same redshift bin for the high 
power AGNs. In Figureoverdensity richness distributions are 
shown for the case of the AGNL sample and for the extractions 
of the corresponding QOL Control sample. The distributions of 
the AGNL sample being sistematically shifted towards higher 
richness values confirm the trend. 

This is another important result of this work: when consider¬ 
ing the environment of radio sources, low power AGNs are found 
preferentially in high density environments. Therefore, higher 
overdensity richness values enhance the probability of a galaxy 
to host a low power radio AGN. This result is not unexpected. 
In fact, it could arise as a direct consequence of different accre¬ 
tion modes for low power and high power AGNs. In particular, 
theory predicts low power radio AGNs to be fuelled by hot gas 
(such as the one found at the center of clusters where cooling 
flows are present), whereas high power AGNs would be fuelled 
by cold gas (see for example Hardcastle et al. 2007, and refer¬ 
ences therein) If this were true, then low power AGNs (fuelled 
by hot gas) would be expected to reside at the centre of rich over¬ 
densities, such as clusters and groups (where hot gas is present). 
On the other hand, high power AGNs would not be found in high 
density environments, both for the lack of cold gas necessary for 
accretion and for the inefficiency of wet mergers (due to the high 
velocity dispersions of galaxies at the centre of clusters and rich 
overdensities), which are the trigger mechanism that has been 
proposed for the ignition of high power radio AGNs. The results 
of this work seem to support this view of AGN and galaxy for¬ 
mation. 


7. Correlation with known ciusters and groups 
cataiogues 

In the previous sections it has been found that radio AGNs live in 
environments which are on average denser than those of normal 
galaxies. It has therefore been explored whether these overden¬ 
sities are related to virialized or bound objects. Two compila¬ 
tions of clusters and groups from the literature have been used: 
that of Finoguenov et al. (2007), found in the COSMOS field 
analysing deep Chandra and XMM observations (see references 
in the paper), and the zCOSMOS group catalogue (Knobel et al. 
2012) For this comparison, only the two most distant redshift 
bins of the AGN sample have been considered (z e [0.7 - I.0[ 
and z G [I.O - 2.0]). At those redshifts, a 0.5 Mpc value cor¬ 
responds approximately to 71 (z = 0.7), 63 (z = 1.0) and 61 
(z = 2.0) arcsec. However, the largest source of uncertainties is 
coming from the large error in the photometric redshifts. 
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To be conservative, all objects that are within a radius of ~ 60 
arcsec and |zlz| < 0.12 for 0.7 < z < 1.0 and ~ 60 arcsec and 
\Az\ < 0.2 for 1.0 < z < 2.0 from the cluster locations reported 
in the two catalogues have been considered, in order to take into 
account the varying error in the photometric redshifts estimate. 

In Table |A.l] the 26 associations with the Finoguenov et al. 
(2007) sample are listed, while in Table |A.2| those with the 
zCOSMOS group catalogue of Knobel et al. (2012) are reported. 
It has been found that ~ 10 - 20% of the AGN sample is in 
groups and clusters. Furthermore, no significant difference be¬ 
tween the radio power distributions of the radio AGNs in clusters 
and groups and of the radio sources of this work is detected. The 
only noticeable source is the one with \og{LiAGHz) - 26.86 in a 
cluster in the most distant redshift bin, which is also the brightest 
object in the AGN sample. 

For what concerns the X-ray luminosities, there is a slight 
tendency of clusters hosting a radio AGN at z € [0.7 - 1.0[ to be 
brighter, with a KS probability of being extracted by the same 
population of ~ 0.09. However the statistics is not sufficient to 
say anything for the highest redshift bin. 

In the case of zCOSMOS groups, the richness distributions 
are significantly different, with groups hosting a radio AGN be¬ 
ing richer (KS probability of ~ 10“^). This holds also when con¬ 
sidering only groups with more than 3 elements. 


8. Mass functions 

It is now investigated whether there is a relative evolution of 
the mass functions of radio AGNs and those of normal galax¬ 
ies. To perform this task, stellar masses as determined in the 
Ultra VISTA Survey (McCracken et al. 2012; Ilbert et al. 2013) 
are used. A limiting magnitude for the UltraVISTA Survey of 
Ks = 24 is assumed in the data. Moreover, the radio power limit 
of\og{Li /^GHz) ^ 24 previously introduced has also been applied 
to the AGN sample (leading to a total of 88 sources), since lack 
of completeness in the radio data could bias the estimate of the 
redshift evolution of the mass functions. 

Stellar mass functions are computed with the non-parametric 
1/Vmax estimator ( Avni & Bahcall 1980, see also Bolzonella 
et al. 2002 and references therein for further information). For 
radio objects, the Vmax is calculated taking into account the dif¬ 
ferent radio and Kg limiting magnitudes of the surveys. In par¬ 
ticular the minimum between the Vmax set by radio and optical 
data has been used. Mass functions have been calculated for all 
UltraVISTA galaxies, a subsample of quiescent ones (selected 
through a cut to \og{SSFRIyr~^) < -II) and the AGN sample. 
The size of the AGN sample is too small to draw conclusions as 
a function of mass, therefore only the integral of the stellar mass 
functions for galaxies more massive than log(M*) = II is con¬ 
sidered (this new cut further reduces the sample to 71 sources). 
These results are shown in Figure[T^ Error bars are computed by 
simply summing all the l/V^^x due to Poisson statistics, 

as explained in Bolzonella et al. (2002) 

Although the integrated mass function for all galaxies varies, 
within a factor ~ 2, between the redshift bins [I.O - 2.0] and 
[0.7 - I.0[ the corresponding variation for the integrated mass 
function of quiescent and radio AGNs is even stronger (a factor 
~ 5 and almost an order of magnitude for the two samples re¬ 
spectively). In these same redshift bins the ratio of radio AGN 
hosting to quiescent galaxies increases from ~ 0.02 to ~ 0.04. 
This increase could confirm that the presence of a radio AGN 
is a phenomenon that evolves with quiescent galaxies, but at a 
different rate. 
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In order to explore if the ratio of radio AGNs to quiescent 
galaxies evolves also as a function of the local density, all sam¬ 
ples have been divided in two, considering objects above and 
below the median of the overdensity richness distributions of 
all galaxies as obtained with our method. With these cuts, the 
numbers of objects in the AGN sample in the high/low density 
regions are, respectively, 13/3, 33/11,17/10 for the low, middle 
and high redshift bins. For the quiescent sample, the statistics is, 
instead, always in the 2000-3000 range in every redshift bin. 

As expected, the objects from the AGN sample are sistemat- 
ically more numerous above the median value than below. This 
could be an indication of the effect previously discovered in a 
more significant way, of radio AGNs being found in higher den¬ 
sity environments. The effect of the environment on the inte¬ 
grated mass functions and on the ratio of AGN hosting to qui¬ 
escent galaxies is showed in Figure [TT] 

The ratio of AGN hosting galaxies to quiescent ones shows 
that the environmental segregation is present. In particular the ra¬ 
tio is larger in high density environments compared to low den¬ 
sity ones at z < 0.7. In the first redshift bin there is a factor 
~ 4 difference between high and low density environments. This 
seems to indicate that the environmental segregation is some¬ 
what more important for radio AGNs. Also, while the evolution 
is similar, both the integrated mass function of the AGN sample 
and the ratio of AGN to quiescent galaxies show an enhanced 
probability of a galaxy to host a radio AGN in high density en¬ 
vironments in the lowest redshift bin. This result is in very good 
agreement with what has been found in previous sections with 
other methods. 

9. Conclusions and summary 

The aim of this work was to investigate the role of the envi¬ 
ronment in determining the presence and the properties of radio 
AGN phenomena. These phenomena are predicted and needed 
by the theory of galaxy formation, therefore understanding the 
way in which environment influences their upcoming could lead 
to a better comprehension of the processes through which galax¬ 
ies are formed. The deep photometric redshift sample extracted 
from the COSMOS survey (Ilbert et al. 2009) has been used, to¬ 
gether with the VLA-COSMOS survey (Schinnerer et al. 2007) 
for the radio data. 

From these surveys, various samples have been extracted, 
with particular focus on AGN samples. Results may be summa¬ 
rized as follows: 

1. Analysing the R sample and the RO Control sample, it is 
found that the environment where radio sources reside is sig¬ 
nificantly denser than the environment around galaxies with¬ 
out radio emission. 

2. These results are not due to radio sources being hosted 
by high mass galaxies (which in turn are known to reside 
in denser environments). In fact, by repeating the environ¬ 
ment estimate on samples designed to have the same stellar 
mass distributions (MR and MO samples) the effect of radio 
sources residing in sistematically and significantly denser en¬ 
vironments still holds. 

3. A sample of radio AGNs has been extracted from the whole 
catalogue of radio sources. It is found that the environment 
around radio AGNs is significantly denser than the environ¬ 
ment around normal galaxies in the same stellar mass and 
Specific Star Formation Rate range (QO sample). This shows 
that the environmental segregation that is found is actually 
due to an enhanced probability of galaxies to host a radio 


AGN in denser environments. This leads to conclude that 
environment plays a role in determining the onset of AGN 
phenomena in galaxies, and that these pehnomena are im¬ 
portant in the process of galaxy formation. 

4. The environmental effect in enhancing the probability of a 
galaxy to host a radio AGN is dominated by low power radio 
AGNs. In fact, the difference in the environment between low 
power radio AGNs and normal galaxies is significant, while 
for high power radio AGNs it is not. This has been found by 
splitting all the previously analysed samples according to a 
cut in their radio luminosity. This effect, of low power radio 
AGNs showing the only signal of environemntal segregation, 
could be related to the thermodynamic properties of the gas 
fuelled to the AGN (see for example Hardcastle et al. 2007) 

5. Analysing the integrated mass functions, calculated for the 
whole population of the UltraVISTA galaxies and sub¬ 
samples of only quiescent and radio AGN hosting ones, hints 
of a redshift evolution can be found. This can be seen both 
in the decrease of the value of the integrated mass functions 
with redshift for the samples of quiescent galaxies and radio 
AGNs, as well as in the same decrease with redshift in their 
ratio. 

6. Exploring the dependence from the environment of the in¬ 
tegrated mass functions, it is found that the fraction of ra¬ 
dio AGNs is always enhanced in high density environments, 
in agreement with the stronger results found using richness 
distributions. Moreover, radio AGNs are also the population 
with the greatest difference in the values of integrated mass 
functions between high and low density environments. 

The picture that emerges from this work is that high density 
environment is indeed related to the presence and the properties 
of radio AGNs. This fact is in good agreement with the current 
theory of galaxy formation, that predicts galaxies in high density 
environments to undergo star formation quenching faster and in 
a more efficient way. Future studies and the increase in the statis¬ 
tics for the mass functions could open the way to a more thor¬ 
ough analysis of the existing correlations. 
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Fig. 5. Galaxy overdensity richness distributions, R and RO sample. 
The top panel refers to 0 < z < 0.7, middle panel to 0.7 < z < 1 
and the bottom panel to 1 < z < 2. The red solid line refers to the 
R sample, the grey lines correspond to the 100 extractions of the RO 
Control sample. The dashed black line is the richness distribution of the 
Control sample extraction corresponding to the median value of the KS 
probability value distribution. 
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Fig. 6. Galaxy overdensity richness distributions, samples MR and MO. 
The top panel refers to 0 < z < 0.7, middle panel to 0.7 < z < 1 
and the bottom panel to 1 < z < 2. The red solid line refers to the 
MR sample, the grey lines correspond to the 100 extractions of the MO 
Control sample. The dashed black line is the richness distribution of the 
Control sample extraction corresponding to the median value of the KS 
probability value distribution. 


Fig. 7. Galaxy overdensity richness distributions, samples AGN and 
QO. The top panel refers to 0 < z < 0.7, middle panel to 0.7 < z < 1 
and the bottom panel to 1 < z < 2. The red solid line refers to the 
AGN sample, the grey lines correspond to the 100 extractions of the 
QO Control sample. The dashed black line is the richness distribution 
of the Control sample extraction corresponding to the median value of 
the KS probability value distribution. 
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Fig. 8. Redshift dependence of radio luminosity. Only the R sample is 
represented. 
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Fig. 9. Galaxy overdensity richness distributions, samples AGNL and 
QOL. The top panel refers to 0 < z < 0.7, the middle panel to 0.7 < 
z < 1 and the bottom panel to 1 < z < 2. The red solid line refers to the 
AGNL sample, the grey lines correspond to the 100 extractions of the 
QOL Control sample. The dashed black line is the richness distribution 
of the Control sample extraction corresponding to the median value of 
the KS probability value distribution. 
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Fig. 10. Redshift evolution of the integrated mass functions. Upper 
panel: value of the integrated mass functions. The curves refer to all Ul¬ 
tra VISTA galaxies (green diamonds), a selection of only the quiescent 
ones (red triangles) and the radio AGN sample (cyan squares). Lower 
panel: ratio of the integrated mass function values for AGN hosting and 
quiescent galaxies. Points have been placed at values of redshift corre¬ 
sponding to the mean redshift of each redshift bin. 



Redshift 


Fig. 11. Environment ejfect on the integrated mass functions. All sam¬ 
ples have been split in two, as stated by the legend, with the terms 
“Over” and “Under” referring to high density environments (over the 
median value) and low density ones (under the median value). Upper 
panel: red and purple triangles refer to quiescent Ultra VISTA galaxies, 
cyan and blue squares to the AGN sample. Lower panel: black and grey 
curves represent the ratio of AGN hosting to quiescent galaxies. Solid 
lines are for high density environments, dashed lines for the low density 
ones. Points have been placed at values of redshift corresponding to the 
mean redshift of each redshift bin. 
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Appendix A: Tables of host X-ray clusters and 
galaxy groups 

In this appendix the complete tables of candidate associations 
between radio AGNs found in this work and known clusters and 
groups catalogues are reported. 
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Table A.l. Candidate associations of radio AGNs and X-ray clusters from the catalogue of Finoguenov et al. (2007) RAcu deCd and Zd are the 
coordinates of the cluster centre. Column 4 is the distance between the cluster centre and the source of the AGN sample considered on the plane 
of the sky, while column 5 is the distance in redshift between the two. Lx([0.1 - 2.A\keV) and \og{Li4GHz) are the X-ray and 1.4 GHz luminosities. 



deCci 

Zd 

Distance (arcsec) 

Az 

Lx([0.1 2A]keV){^Z^) 

1 / Li^ACHz 

0.7 < z < 1 

150.41386 

1.84759 

0.969 

40.176 

0.015 

2.9883 X 10« 

25.68 

150.27736 

2.05303 

0.908 

41.112 

-0.0521 

9.6268 X 10^2 

23.76 

150.02382 

2.20323 

0.942 

7.1856 

0.0988 

1.9325 X 10^3 

23.35 

149.64966 

2.20925 

0.954 

0.120744 

-0.0165 

1.7111 X 10^3 

24.34 

150.21454 

2.28010 

0.881 

26.3124 

-0.0419 

1.3734 X 10^3 

24.59 

149.95262 

2.34188 

0.942 

0.0144 

0.0156 

1.5213 X 10^3 

23.47 

149.95262 

2.34188 

0.942 

41.76 

-0.0002 

1.5213 X 10^3 

24.16 

149.92926 

2.40902 

0.874 

0.0072 

-0.0044 

8.4425 X 10^2 

23.58 

149.63988 

2.34912 

0.951 

41.112 

-0.0369 

3.2975 X 10^2 

24.06 

149.62355 

2.39918 

0.845 

0.108072 

0.0051 

1.9522 X 10^2 

24.31 

150.15298 

2.39447 

0.899 

26.01 

0.0193 

9.4442 X 10^2 

23.64 

149.66927 

2.47365 

0.957 

11.4408 

-0.0044 

1.8343 X 10^2 

23.92 

149.56174 

2.42195 

0.846 

0.0144 

0.0332 

1.7645 X 10^2 

24.20 

150.00713 

2.45343 

0.731 

0.1476 

0.0229 

5.2795 X 10^2 

24.15 

149.92343 

2.52499 

0.729 

30.8016 

0.0016 

1.1257 X 10^^ 

23.74 

150.10533 

2.72392 

0.727 

2.39364 

-0.0023 

6.9659 X 10^2 

23.89 

149.91772 

2.70088 

0.889 

3.00132 

-0.0355 

2.1160 X 10^2 

24.28 

150.58397 

2.32155 

0.720 

46.116 

0.0182 

4.2323 X 10^2 

23.39 

150.05057 

2.13923 

0.959 

1.09008 

0.1189 

5.7184 X 10^2 

24.19 

1 < z < 2 

150.76245 

1.79362 

1.258 

10.4004 

-0.0051 

5.5689 X 10« 

26.43 

149.51855 

2.09959 

1.382 

48.096 

0.1839 

7.8229 X 10^2 

24.00 

150.59309 

2.53890 

1.045 

7.056 

-0.1063 

3.1937 X 10^2 

23.63 

149.59763 

2.44004 

1.168 

10.4076 

-0.023 

2.6946 X 10^2 

26.86 

150.57024 

2.49864 

1.146 

17.586 

-0.0083 

2.5386 X 10^2 

24.13 

150.12646 

1.99926 

1.019 

13.644 

-0.1744 

1.2897 X 10^2 

23.65 
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Table A.2. Candidate association between radio AGN and groups from the catalogue of Knobel et al. (2012) RAgr, decgr and Zgr are the coordinates 
of the group centre. N spec IN photo is the number of sources with spectroscopic and photometric redshift respectively, column 5 and 6 are the distances 
from the group centre to the position of the radio AGN source on the plane of the sky and along the redshift direction respectively. \og{Li_^GHz) is 
the 1.4 GHz luminosity. 


RAgr 

dec gy 

H spec / H photo 

Zgr 

Distance(arcsec) 

Az 

1 / LrtGHz X 

0.7 < z < 1 

149.920567 

2.521800 

12/82 

0.7297 

17.6328 

0.0023 

23.74 

150.005475 

2.451841 

12/49 

0.7311 

8.1792 

0.023 

24.15 

150.212666 

2.281762 

7/14 

0.8812 

18.7812 

-0.0417 

24.59 

149.552005 

2.423054 

5/25 

0.8446 

35.2692 

0.0318 

24.2 

149.914194 

2.694681 

5/54 

0.8898 

55.008 

-0.0669 

24.05 

149.914194 

2.694681 

5/54 

0.8898 

28.5552 

-0.0347 

24.28 

150.101070 

2.268664 

3/4 

0.6855 

47.16 

-0.0547 

23.36 

150.094158 

2.063360 

311 

0.725 

53.064 

0.01 

23.61 

150.580055 

2.328653 

3/6 

0.7274 

28.3644 

0.0256 

23.39 

150.023066 

2.516794 

3/2 

0.7473 

31.554 

0.0314 

24.92 

150.508547 

2.653296 

311 

0.8102 

34.056 

0.0103 

25.02 

149.569460 

2.418630 

3/107 

0.8529 

30.258 

0.0401 

24.2 

150.515963 

2.005564 

3/11 

0.8797 

15.822 

0.0089 

24.29 

149.788439 

2.757904 

3/10 

0.9093 

24.6348 

-0.0272 

24.16 

150.475201 

1.626689 

3/1 

0.9713 

32.6952 

-0.0236 

23.64 

150.098295 

2.056180 

2/3 

0.6375 

47.34 

-0.0775 

23.61 

150.017350 

2.442392 

2/5 

0.6681 

54.288 

-0.04 

24.15 

149.725800 

2.770682 

2/7 

0.7032 

42.48 

-0.0508 

23.98 

149.838840 

1.683059 

2/7 

0.747 

0.126828 

0.0067 

23.3 

149.731350 

2.759319 

2/6 

0.7655 

18.3024 

0.0115 

23.98 

150.231175 

2.072644 

2/1 

0.7995 

40.212 

-0.0594 

23.61 

150.465750 

2.423217 

2/4 

0.8266 

9.9792 

-0.0006 

24.4 

150.410760 

1.801465 

2/5 

0.8456 

13.6152 

0.0372 

24.77 

149.985860 

2.325898 

2/3 

0.849 

32.7384 

-0.0764 

23.94 

150.091666 

2.599807 

2/4 

0.8929 

20.2572 

-0.0961 

24.07 

1 < z < 2 

150.216932 

2.273971 

2/5 

0.893 

44.424 

-0.0299 

24.59 

150.445160 

1.845406 

2/8 

0.8943 

20.6964 

0.0232 

24.6 

150.023318 

2.205914 

2/12 

0.9409 

5.7276 

0.0977 

23.35 

149.963878 

2.363113 

2/10 

0.9454 

50.04 

0.0032 

24.16 

149.649660 

2.209250 

2/7 

0.9539 

0.120744 

-0.0166 

24.34 

149.843245 

2.573333 

2/8 

0.9657 

42.408 

-0.0285 

25.08 

149.633185 

2.457927 

3/2 

1.1688 

40.716 

-0.0832 

25.06 

150.128820 

1.922544 

2/1 

1.0083 

53.172 

-0.0832 

23.64 










